Observational tests of a two parameter power-law class modified gravity in Palatini 

formalism 
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CONTEXT: In this work we propose a modified gravity action f(R) = (R n - R%) 1,n with two 
free parameters of n and Rq and derive the dynamics of a universe for this action in the Palatini 
formalism. 

AIM: We do a cosmological comparison of this model with observed data to find the best parame- 
ters of a model in a flat universe. 

METHOD: To constrain the free parameters of model we use SNIa type la data in two sets of gold 
and union samples, CMB-shift parameter, baryon acoustic oscillation, gas mass fraction in cluster 
of galaxies, and large-scale structure data. 

RESULT: The best fit from the observational data results in the parameters of model in the range 
of 7i = 0.98±g;g| and Q, M = 0.25+J};} with one sigma level of confidence where a standard ACDM 
universe resides in this range of solution. 

PACS numbers: 04.50. +h, 95.36,+x, 98.80.-k 
I. INTRODUCTION 

Recent observation of CMB+SNIa reveals that the 
Universe is under positive acceleration, in contrast to our 
expectations from the behavior of ordinary matter. One 
of the possible solutions is assuming a cosmological con- 
stant to provide a late time acceleration to the universe 
[1]. Although ACDM is the easiest model that fits well 
with the current observational data, it suffers from the 
fine-tuning and coincidence problems, which motivates to 
introduce alternative theories such as dark energy mod- 
els. The other possibility is the modification of grav- 
ity law in such a way that behaves as standard general 
relativity in a strong gravitational regime, but repulses 
particles in the cosmological scales [2]. Each proposed 
model must challenge two different observational criteria 
of (a) cosmological tests and (b) local gravity tests in 
solar system scales. 

In this work we examine a f(R) = (R n — R$ ) 1 ' n gravity 
model with the cosmological observational data to fix the 
two parameters of the model, i?o and n. We use Super- 
nova Type la, CMB shift parameter, baryonic acoustic 
oscillation, gas mass fraction in cluster of galaxies and 
large scale structures to find the best parameter of the 
model. The organization of paper is as follows: In Sec. 
II, we introduce the field equation in f(R) modified grav- 
ity model in Palatini formalism. Using FRW metric we 
derive a modified Friedman equation. In Sec. Ill we ap- 
ply our proposed action to derive the field equations. In 
Sec. IV we study geometrical and dynamical behavior of 
universe in this model. In Sec. V using the geometri- 
cal observations as SNIa, CMB-shift parameter, baryon 
acoustic oscillation, gas mass fraction of clusters of galax- 



ies, we put constraints on the parameters of the model. 
Finally, we do comparison of models with the large scale 
structures' data in Sec. V. The conclusions arc given in 
Sec. VI. We show that the best parameter of this model 
privilege a ACDM universe. 



II. MODIFIED GRAVITY IN PALATINI 
FORMALISM 



For an arbitrary action of the gravity as a function of 
Ricci scalar f(R), there are two main approaches to ex- 
tract the field equations. The first one is the so-called 
metric formalism, which is obtained by the variation of 
action with respect to the metric. In this formalism in 
contrast to the Einstcin-Hilbert action, the field equation 
is a fourth order nonlinear differential equation. In the 
second approach so-called Palatini formalism, the con- 
nection and metric are considered independent fields and 
variation of action with respect to these fields results in 
a set of second order differential equations. In what fol- 
lows we will work in Palatini formalism. Let us take the 
general form of action in the Palatini formalism as 

S[f;g,t^ m } = ^ J d 4 x^f(R) + S m [9^,*m}, (1) 

where k = 8ttG and S , m [9/ii/, 4"m] is the matter action 
that depends on metric g^ v and the matter fields ^> m . 
R = R(g,T) = g^ v R^ v {t) is the generalized Ricci scalar 
and R^ v is the Ricci tensor, made of affine connection. 
Varying action with respect to the metric results in 

1 



f'(R)R flv (T) - -f(R)g, 



kT,, 



(2) 



where prime is the differential with respect to the Ricci 
scalar and is the energy-momentum tensor 



T — 



-2 SS m 



(3) 



On the other hand varying the action with respect to the 
connection results in 



(4) 



where V is the covariant derivative defined from parallel 
transformation and depends on affinc connection. From 
Eq.(4), we can define a new metric of = f{R)g^.i, 
conformally related to the physical metric where the con- 
nection is the Christoffcl symbol of this new metric. We 
take a flat FRW metric (namely K = 0) for the universe 



-dt 2 



a{t) 2 8ijdx 1 'dx\ 



(5) 



and assume that universe is filled with a perfect fluid with 
the energy-momentum tensor of T" = diag(—p,p,p,p). 
Using the metric and energy momentum tensor in Eq.(2) 
we obtain the generalized FRW equations. It should be 
noted that the conservation law of energy-momentum 
tensor, T^ v = is defined according to the covariant 
derivative with respect to the metric to guarantee the 
motion of particles on geodesies [3]. A combination of 
GS and G) results in 



1 k{ P + 3p) If 
6 /' ' 



2f" 6 /' 
On the other hand, the trace of Eq. (2) gives, 
Rf'(R) - 2f(R) = K T, 



(G) 



(7) 



where T = g^T^ = —p + 3p. The time derivative of 
this equation results in R in terms of the time derivative 
of density and pressure. Using the equation of state of 
cosmic fluid p = p(p) and continuity equation, the time 
derivative of Ricci is obtained as 



R = 3kH 



(l-3dp/dp)(p + p) 
Rf" - f'(R) 



(8) 



To obtain a generalized first FRW equation, we start with 
Eq. (7) and obtain the density of matter in terms of the 
Ricci scalar as 



up 



2/ - Rf 
1 -3w ' 



(9) 



where w = p/p. Substituting Eq. (9) in (6) and using Eq. 
(8) to change d/dt = Rd/dR, we obtain the dynamics of 
the universe in terms of the Ricci scalar as 



H 2 



1 



3(l + u)f-(l + 3u)Rf 



6(1 -3u)f 
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+ f(l + w)^ 



f"{2f-Rf>) 
'(ft/"-/') 



(10) 



On the other hand using Eq. (7) and the continuity equa- 
tion, the scale factor can be obtained in terms of the Ricci 
scalar 



Kpo(l — 3cj) 



(2/ - Rf) 



(11) 



where po is the energy density at the present time and 
ao, the scale factor at the present time, is set to 1. Now 
for a generic modified action, eliminating the Ricci scalar 
in favor of the scale factor between Eqs. (10) and (11) we 
can obtain the dynamics of universe [i.e. H = H(a)]. For 
the simple case of matter dominant epoch uj = 0, these 
equations reduce to 



H 1 



1 



3/ - Rf 



6f 



3 f"(2f-Rf) 
2 /'(H/"-/') 



(12) 



and 



(2/ — Rf) 

Kp 



(13) 



III. F(R) = (R N - R^) 1/N GRAVITY 

Here in this section we propose a modified gravity ac- 
tion of f(R) = (R n - i? ™) 1/ " with the two free parame- 
ters where Rq > and n > 0. This action is a generalized 
form of n = 2 that has been discussed in [4,5]. This ac- 
tion has a minimum vacuum in an empty universe and a 
flat Minkowski space is not achievable in this action. This 
behavior causes an accelerating expansion of the universe 
for a low density universe. The minimum curvature from 
the vacuum solution in Eq. (7) is: 



R v — 2 1 /"i?o- 



(14) 



On the other hand, for the strong gravitational regimes 
the action reduces to the Einstein-Hilbert action. To 
have the asymptotic behavior of action for these two ex- 
treme cases we do a Taylor expansion of action around 
R v in an almost empty universe and R v / R — ► at strong 
gravitational regimes. For the weak field, the expansion 
of action results in 



f(R) = R v (\)^ + (\^- 1 (R-R v ) 
+ { 1 - ) i^l(R-R v r + ... 



(15) 



where ignoring higher order terms, we can rewrite this 
equation as 



f(R) = R-A(R v ,n). 



(16) 



Here A(R v ,n) is an effective cosmological constant de- 
pends on the curvature in a vacuum and the exponent 



of action. On the other hand we expand the action in a 
strong gravitational regime (e.g. R^> R v ). In this case 
the action can be written as follows: 



m— 1 



* 1 1 D 

f(R) = R+J2 — U(-- k)(-ir(^-r n R. (l?) 

^-^ m! - L - L n ft 

m=l fe=0 

Ignoring the higher orders in a strong gravitational field, 
this action reduce to the Einstcin-Hilbert action. So our 
chosen actions in these two extreme regimes vary from 
the Einstein-Hilbert to the Einstcin-Hilbert plus cosmo- 
logical constant. 

Let us study the solutions of modified gravity in three 
common cases of a pointlike source in vacuum space, a 
universe in radiation, and matter-dominant epochs. For 
a pointlikc source in an empty space, letting T^ v = 
outside the star, wc will have a constant Ricci scalar. 
This means that we will have constant R v , f(R v ) and 
f'(Rv) for all the space. With this condition we write 
the field equation as follows: 



r f'(Rv) >9flv ' 



(18) 



where the coefficient of metric at the right-hand side 
of the equation plays the role of effective cosmological 
constant. Substituting the corresponding value for the 
vacuum from Eq. (14), the effective cosmological con- 
stant obtain A e // = 2~~ 2 Rq. So the solution of the field 
equation in the spherically symmetric space results in a 
Schwarzschild-de-Sitter space. The value of Ro will be 
fixed in the next section from the cosmological observa- 
tions. 

In the radiation dominated era wc have p = p/3. With 
this equation of state, the trace of the energy-momentum 
tensor is zero and it resembles a vacuum solution where 
the Ricci scalar is constant and equal to 2 1//n i?o- We 
substitute f(R) and its derivatives in Eq.(6) to have the 
dynamics of the Hubble parameter as a function of den- 
sity of the universe 



E z = 



-(2np + R a 



(19) 



Analysis for n = 2 shows that Ro is in the order of Hq 
[4,5]. So for the radiation-dominant epoch, we neglect 
Ro in comparison with the density of the universe. Using 
the continuity equation provides p cx a~ 4 , then the scale 
factor changes with time as a a t 1 ! 2 . This result shows 
no dynamical deviation from the standard cosmology at 
the early universe. 

For the matter-dominant epoch, we calculate the dy- 
namics of the universe for simplicity in terms of a new 
variable, X = R/H 2 . The action can be written in this 
new form as 



f(R) = H Q 2 F(X) 
F(X) = {X n - X n ) 1/n 



(20) 
(21) 



where X = Ro/Hg and H = lOO/i Km/s/Mpc. The 
relation between the derivatives with respect to R and 
new variable X is related as 



f(R) = F'(X), 
F"(X) 



f"(R) 



Hi 



(22) 
(23) 



where the derivatives in the left-hand side of equations 
are with respect to the Ricci scalar, but in the right-hand 
side they are in terms of X , (i.e. ' = t^f)- 

We rewrite Eq. (12) with the new dimensionlcss pa- 
rameter X: 



U 2 (X) 



1 



3F - XF 1 



6F' 



(1 



3 F"(2f-AT) „ : 
2 F'(XF"-F'j I 



(24) 



where Ti{X) = H/Hq is the normalized Hubble parame- 
ter to its current value. Using the conventional definition 
of Q m at the present time as VL„ 
(9) we obtain 



np^/SH'l and Eq. 



n m (x) = 



2F - XF' 



(25) 



where f2 m (A) = Cl m a 3 . We obtain the relation between 
the scale factor and dimensionlcss parameter X as: 



2F - XF' 

3 ^ m 



-1/3 



(26) 



In order to have positive scale factor, X should change 
in the range of X > 2~Xq, where the minimum value for 
Xq is in agreement with the vacuum solution of the Ricci 
scalar. 

An important point worth mentioning here is that the 
two free parameters Xq and n, appeared in the dynamics 
of the universe can be replaced with more relevant ones. 
One of them is f2 m presented in Eq.(25), where depends 
directly to n, X and X p (p stands for the present time). 
X p can be eliminated using Eq. (24), letting TC(X p ) = 1 
results in a relation between Xp and Xq and n. The 
second parameter we will use instead of Xq is X p . 



IV. GEOMETRICAL PARAMETERS IN F(R) 
GRAVITY 



The cosmological observations are mainly dependent 
on background spatial curvature and four dimensional 
space-time curvature of the Universe. In this section we 
introduce the geometrical parameters in modified gravity 
to use it in observational tests of model. 



A. comoving distance 

The radial comoving distance is one of the basic pa- 
rameters in cosmology. For an object with a rcdshift of 
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FIG. 1. The dependence of comoving distance as a function 
of redshift for the case of n — 1 and various Xq. Increasing 
Xo makes larger comoving distance for a given redshift. 
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FIG. 2. The dependence of comoving distance as a func- 
tion of redshift for the case of Xo = 4.2 and various n. In 
increasing the exponent results in smaller comoving distance 
for a given redshift. 



z, using the null geodesies in the FRW metric, the co- 
moving distance in terms of X is obtained by 



r(z;n,X ) = c 



dz' 

H( z y 



(27) 



cm 



F' - XF" dX 
3 4 /3(tt m )V3 J Xp (2F-XF')iH(X)' 

(28) 

where the dimensionless parameter X relates to the red- 
shift from equation (11) as: 

z = - X o) " _1 (* n - 2X o )] 1/3 " 1 ( 29 ) 

OS L m 

Knowing the parameters of the action n and Xq, we 
can calculate the Hubble parameter at a give X by Eq. 
(24) substituting it in (28) we obtain a comoving distance 
by numerical integration. Figures (1) and (2) show co- 
moving distance as a function of redshift in the unit of 
cHq 1 for various values of parameters of the model. In 
Fig. (1) we fix n = 1 which is equivalent to the ACDM 
universe and let Xo vary. It seems that Xq plays the role 
of effective cosmological constant. Increasing this term 
makes a larger comoving distance for a given redshift. In 
Fig. (2) we keep X = 4.2 and let n change. Increasing 
the exponent results in a smaller comoving distance for 
a given redshift. 



B. Angular diameter distance and Alcock-Paczynski 

test 

The apparent angular size of an object located at the 
cosmological distance is another important parameter 
that can be affected by the cosmological model. An ob- 
ject at the redshift of z and the perpendicular size of D± 
is seen by the angular size of 



A6 



El 

dA ' 



(30) 



where dA = r{z\n,Xo)/{\ + z) is the angular diameter 
distance. Now imagine this structure has the size of Dii 
along our line of sight. Then the light arriving at us from 
the back and front of this structure will not have the same 
redshift. The difference in the redshifts of the two sides 
of the structure can be obtained by the delay in received 
light to the observer with At(z) = Di\/c. Writing At in 
terms of Aa as At = H~ 1 (z)Aa/a we again change Aa 
in terms of Az as Az/(1 + z) = —Aa/a. The result is 
writing the width of the structure in redshift space along 
our line of sight in terms of physical size as 



Az = -DuH(z)(l + z). 
c 11 



(31) 



Now the width of the structure in the redshift space to 
the apparent angular size of structure obtain as 



Az = (l + z)H(z)d A ( D\\ . 
A6 c { Dj' 



(32) 
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FIG. 3. The dependence of Az/ AO as a function of redshift 
for the case of n = 1 and various Xq. Increasing Xo causes 
increasing the apparent size of cosmological objects. 



For the spherical structures for instance taking into ac- 
count the neutral hydrogen clouds at z < 6 with the 
spherical symmetric shape, Eq. (32) is written as: 

Az _ H(z;n,X Q )r(z;n,X ) 
A9 ~ c ' 

This relation is the so-called Alcock-Paczynski test. The 
advantage of the Alcock-Paczynski test is that this re- 
lation is independent of the Hubble parameter at the 
present time and of the existence of the dust in the in- 
tergalactic medium. In this method, instead of using a 
standard candle, we will use a standard ruler such as the 
baryonic acoustic oscillation. 

Figures (3) and (4) show a dependence of Az/ AO as 
a function of redshift normalized to the corresponding 
value in a ACDM universe. In Fig. (3) the relative size 
of the structure in redshift space to the observed angular 
size is compared to that in a ACDM universe for a fixed 
value of n = 1. In Fig. (4) we fixed Xq = 4.2 and change 
the exponent n. 

C. Comoving Volume Element 

The comoving volume element is another geometrical 
parameter that is used in number-count tests such as 
lensed quasars, galaxies, or clusters of galaxies. The co- 
moving volume element in terms of comoving distance 
and Hubble parameter is given by 
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FIG. 4. The dependence of Az/A9 as a function of red- 
shift for the case of Xo = 4.2 and various n. Increasing the 
exponent results in decreasing the apparent cosmological size 
of objects. 



TABLE I. Different priors on the parameter space, used in 
the likelihood analysis. 



Parameter 


Prior 




K 


0.00 


Fixed 


Q, b h z 


0.020 ± 0.005 


Top hat (BBN) [6] 


h 




Free [7,8] 


w 





Fixed 



Figures (5) and (6) show the dependence of comoving 
volume clement as a function of redshift. Figure (5) rep- 
resents the dependence of comoving volume for a fixed 
n = 1 and various Xq. Increasing Xq causes larger co- 
moving volume element. In Fig. (6) we plot the volume 
for fixed Xo = 4.2 changing n. Increasing the exponent 
index makes the comoving volume element smaller. 

V. OBSERVATIONAL CONSTRAINTS: 
BACKGROUND EVOLUTION 

In this section we compare the observed data with that 
from the dynamics of the background from the model. 
We use Supernova Type la data, CMB-shift parameter, 
baryonic acoustic oscillation (BAO) and the gas mass 
fraction of cluster of galaxies to constrain the parameters 
of the model. 
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FIG. 5. The dependence of comoving element as a function 
of redshift for the case of n — 1 and various Xq. Increasing 
Xq increases the comoving volume element. 
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FIG. 6. The dependence of comoving element as a function 
of redshift for the case of Xq = 4.2 and various n. Increasing 
the exponent causes decreasing the comoving volume element. 



A. Supernova Type la 

The Supernova Type la experiments provided the main 
evidence for the existence of dark energy. Since 1995 two 
teams of High-Z Supernova Search and the Supernova 
Cosmology Project have discovered several type la su- 
pernovas at the high redshifts [9,10]. They showed that 
to interpret the faintness of high redshift supernovas in a 
flat universe one has to consider an accelerating universe 
at the present time. 

In this work we take two sets of SNIa data. The first 
one is the gold sample which has a 157 supernova [11] and 
the second set is a combined data set of a 192 supernova 
[12]. The distance modulus for supernovas is calculated 
by 

/i = m — M = 5 log Dl(z; Xq, n) 

c 



where 



D L (z;X Q ,n) = (l + z)H 



dz' 
~H&~V 



(36) 



and Dl can be written in terms of new parameter 
AT, which appeared in the redefinition of the modified 
gravity action as 

_ 1 {2F - XF')i f x F'-XF" dX 
L 3 ( 3 fi m )l J Xp (2F-XF')iH(Xy [ > 

For simplicity in calculation, we define 



(38) 



which is a function of the Hubble constant at the present 
time. We write the distance modulus as 



H = 5logD L {z;X ,n) + M. 



(39) 



In the next step we use \ 2 fitting to constrain the pa- 
rameters of the model. 



X 2 (M,X 0l n)=J2 



[ptobs(zi) - /j 4h (z i ;M,X ,n,)}' i 



(40) 



where <7j is the uncertainty in the distance modulus. To 
constrain the parameters of the model, we use the likeli- 
hood statistical analysis 



C(M,X ) =AAe^ 2 (A^o)/ 2i 



(41) 



where Af is a normalization factor. The parameter M is 
a nuisance parameter and should be marginalized (inte- 
grated out) leading to a new \ 2 defined as: 



X z = -2 In / e~ x ~ /2 dM. 

J — OO 

Using Eqs. (40) and (42), we find 



(42) 



X 2 (X ) = X 2 (M = 0, Xo) ~ + ln(C/2vr), (43) 



C 



where 




\p*obs(zi) - Vth{zi;X ,M = 0)] 



and 



(44) 



(45) 



Equivalent to marginalization is the minimization of x 2 
with respect to M. One can show that \ 2 can be ex- 
panded in terms of M : 

X 2 sm,(Xo) = X 2 (M = 0, X ) - 2MB + M 2 C\ (46) 

which has a minimum value for M = B/C and results 



xlm a (Xo)=X 2 (M = 0,X ) 



C 



(47) 



Using equation (47) we can find the best fit values of 
model parameters, minimizing XsNia(^o)- Figures (7) 
and (8) represent the best fit to the Supernova Type la 
new gold sample and union sample respectively. The best 
fit values for the free parameter of the model for two cases 
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the new gold sample and n 



l-63i ' 92 , X 



7+0.64 

64 

1+1.32 

2.86 



0.67™ for 



6.09; 



and VL m = O.47+0 47 for mixed Gold-SNLS data. Figures 
(9) and (10) represent the likelihood functions in terms 
of n and Xq . 



B. CMBR Shift parameter 

Another dynamical parameter that is used in recent 
cosmological tests is the CMB shift parameter. Before 
the last scattering epoch, the baryons and photons were 



tightly coupled through the electromagnetic interaction. 
This coupled fluid was under the influence of two ma- 
jor forces of (a) the gravitational pull of matter and (b) 
the out leading pressure of photons. The finger print of 
this competition leads to the familiar spectrum of peaks 
and troughs on the CMB map. Here the main peak is 
the so-called acoustic peak. The odd peaks of the CMB 
anisotropy spectrum correspond to the maximum com- 
pression of the fluid, the even ones to the rarefaction 
[13]. 

In an idealized model of the fluid, there is an analytic 
relation for the location of the m-th peak: l m w ml a 
[14,15], where I a is the acoustic scale which may be cal- 
culated analytically and depends on both pre- and post- 
recombination physics as well as the geometry of the uni- 
verse. The acoustic scale corresponds to the Jeans length 
of photon-baryon structures at the last scattering surface 
some ~ 379 Kyr after the big bang [16]. The appar- 
ent angular size of the acoustic peak can be obtained by 
dividing the comoving size of the sound horizon at the 
decoupling epoch r s (z c i ec ) by the comoving distance of 
observer to the last scattering surface r(zd ec ): 



TV 

Ta 



(Zdec) 



r(zdec) 



(48) 



The nominator of Eq. (48) corresponds to the distance 
that the perturbation of pressure can travel from the big 
bang to up to the last scattering surface, which is defined 
as the integral below: 



° v s (z')dz' 
ec H(z')/H 



(49) 



where v s (z) 2 = 3 + 9/4 x pb{z) / p ra d{z) is the sound 
velocity in the unit of speed of light from the big bang 
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FIG. 9. Marginalized likelihood functions of f(R) modified 
gravity model free parameter, n. The solid and dash lines cor- 
respond to the likelihood function of fitting the model with 
SNIa data new gold sample and union data set, respectively. 
The intersections of the curves with the horizontal solid and 
dashed lines give the bounds with la and 2a level of confi- 
dence, respectively. 
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FIG. 10. Marginalized likelihood functions of f(R) mod- 
ified gravity model free parameter, Xq. The solid and dash 
lines correspond to the likelihood function of fitting the model 
with SNIa data new gold sample and union data set, respec- 
tively. The intersections of the curves with the horizontal 
solid and dashed lines give the bounds with la and 2a level 
of confidence, respectively. 



up to the last scattering surface [14,17] Zd ec is the redshift 
of the last scattering surface. 

Changing the parameters of the model can change the 
size of the apparent acoustic peak and subsequently the 
position of I a = "k/6a i n the power spectrum of tempera- 
ture fluctuations on CMB. The simple relation l m ss ml a 
however does not hold very well for the first peak al- 
though it is better for higher peaks [2]. Driving effects 
from the decay of the gravitational potential as well as 
contributions from the Doppler shift of the oscillating 
fluid introduce a shift in the spectrum. A good param- 
eterizations for the location of the peaks and troughs is 
given by 



lm = lA{m~ (j) m ) 



(50) 



where <p m is a phase shift determined predominantly by 
prerecombination physics, and is independent of the ge- 
ometry of the Universe. Instead of the peak locations of 
the power spectrum of CMB, one can use another model- 
independent parameter, which is the so-called shift pa- 
rameter 72, as 



1/2 



72 



Yj^sinn k (uj k r) 



(51) 



where sinnk{x) — sin(a;), x, sinh(x) for k = —1,0,1. For 
the case of a flat universe, which is our concern, the shift 
parameter reduces to the simpler formula of 



72 



dz 



(52) 



Now we change the variable from the redshift to X and 
rewrite the above expression in terms of dimcnsionlcss 



parameter X and take the integral from the value of X 
at the present time as the lower limit of the integral and 
the value of X at the decoupling time as the upper limit: 



72 



- - r; 



X, 



pi _ XF" 



dX 



3* Jx p (2F-XF')iH(X) 



(53) 



The observed result of the CMB experiment is 72 = 
1.716 ± 0.062 [16]. It is worthwhile to mention that the 
dimensionless parameter 72 is independent of the Hub- 
ble constant. We compare the observed shift parameter 
with that of the model using the likelihood analyzing, 
minimizing x 2 defined as 

Xcmb = [7 "° b % W (54) 

°CMB 



C. Baryon Acoustic Oscillations 

Another geometrical cosmological probe which deter- 
mines the distance- redshift relation is BAO. The physics 
governing the production of BAO is well understood. 
Acoustic peaks occurred because cosmological perturba- 
tions excite sound waves in initial relativistic plasma in 
the early epoch of the Universe. Dark matter perturba- 
tions grows in place while the baryonic matter pertur- 
bations were carried out in an expanding spherical wave 
because of their interaction with photons. At the recom- 
bination epoch, when the photons started to decouple 
from the baryonic matter, the shell of the baryonic mat- 
ter perturbation sphere was nearly 150 Mpc in the co- 
moving frame. From the linear structure formation the- 
ories, this scale should not be changed until the present 



time. The structure formation theory predicts that this 
150Mpc imprint of baryonic matter remains in the cor- 
relation function of the density contrast and can be seen 
in the large-scale surveys. 

By knowing the size of acoustic oscillation, one can 
measure the angular distance to this structure. The 
large-scale correlation function measured from 46748 lu- 
minous red galaxies spectroscopic sample of SDSS in- 
clude a clear peak at 100 Mpcft. -1 [18]. The correspond- 
ing comoving scale of the sound horizon shell is about 
150Mpc in radius. A dimcnsionless and Hq independent 
parameter for constraining the cosmological models has 
been proposed in literatures [18] as follows: 




HaD 2 L (z s & ss ; X ) 



H(z sdss ; X )z 2 dss (l + z sdss y 



1/3 



(55) 



(56) 



We rewrite the above dimcnsionless quantity in terms of 
modified gravity model parameters as 

'(3Q TO )- 



1 


p d3B dz 


_ ^sdss J 


o H(X)\ 



3z s 



A' , 



(F' - XF")dX 
Xp H(X)(2F - XF') 2/3 



(57) 



Now we can put a constraint on the f(R) modified grav- 
ity model using the value of A = 0.469 ± 0.017 from 
luminous red galaxies observation at zsdss = 0.35 [18]. 
It is worthwhile to mention that the procedure above 
presented in literature is well proposed for dark energy 
models in which the Q m has the same definition in stan- 
dard cosmology. In contrast to the dark energy models in 
gravity theories in the Palatini formalism, fl m is a con- 
ventional dimcnsionless parameter and does not have the 
same role as in dark Energy models, considering the well- 
known fact that Q m = 1 does not correspond to a flat 
universe as in standard FRW equations. Consequently in 
order to not include the weak model dependence of the 
dimensionless parameter, A, we will use another similar 
approach, proposed by Percival ct al. [19]. This method 
constrain general cosmological models by using BAO dis- 
tance measurement from galaxy samples covering differ- 
ent redshift ranges. Measuring the distance redshift rela- 
tion at two rcdshifts of z = 0.2 and z = 0.35 for clustering 
of SDSS luminous red galaxies enables us to define a new 
dimcnsionless parameter as 



B = 



D v (z = 0.35) 
D v (z = 0.20) 



where Dy is given by 



D v = 



(1 + zfd\cz 



1 V3 



H(z) 



(58) 



(59) 



and dA is the angular diameter distance. The observa- 
tional values for two different redshifts are reported in 
[19] with la error: 



D v (z = 0.20) 
D v (z = 0.35) 



= 0.1980 ± 0.0058 
= 0.1094 ±0.0033 



(60) 

(61) 
(62) 



where r s is the comoving sound horizon scale at the re- 
combination epoch. Considering that BAO measurements 
have the same measured scale at all redshifts then we 
have a numerical value for B as 



B 



Dv{z = 0.35) 
D v (z = 0.2) 



1.812 ±0.060. 



(63) 



Now we convert the comoving angular diameter distance 
to luminosity distance Dl and calculate Dl in terms of 
the dimcnsionless Hubble parameter in modified gravity 



B 



H{z = 0.20)D L 2 (z = 0.35)0.35(1 + 0.2) ; 



,1/3 



H(z = 0.35)D L 2 (z = 0.20)0.2(1 + 0.35) 2 



(64) 



We use x 2 as one more fitting parameter with the ob- 
served value of B = 1.812 ± 0.060. This observation per- 
mits us to add one more term to \ 2 from that of SNIa 
and CMB-shift parameter by minimizing 



2 

Xbao 



{B obs - B th f 



(65) 



BAO 



This is the third geometrical parameter we will use to 
constrain the model. 



D. Gas mass fraction of cluster of galaxies 

Measurement of the ratio of X-ray emitting gas to the 
total mass in galaxy clusters (f ga s) a lso is an indication 
of the acceleration of the Universe. This method can be 
used as another cosmological test to constraint the pa- 
rameters of the model. Galaxy clusters are the largest 
objects in the Universe; the gas fraction in them is pre- 
sumed to be constant and nearly equal to the baryon 
fraction in the Universe. Sasaki (1996) and Pen (1997) 
described how measurements of the apparent dependence 
of the baryonic mass fraction could also, in principle, be 
used to constrain the geometry and matter content of a 
universe [20,21]. The geometrical constraint arises from 
the dependancc of the measured baryonic mass fraction 
value on the assumed angular diameter distance to the 
clusters [22]. The baryonic mass content of galaxy clus- 
ters is dominated by the X-ray emitting intcrclustcr gas, 



the mass of which exceeds the mass of optically luminous 
material by a factor of 6 [23,24]. Let us define f gas as 



fgas 



M,, 



(66) 



In the second step we want to replace the mass of gas by 
the baryonic mass considering that : 



M b = (1 + P)M ga 



gas i 



(67) 



where from the observations we know [3 = O.Wh 1 / 2 [23]. 
On the other hand we assume that we are observing rich 
cluster of galaxies where the fraction of baryonic mass to 
the total mass has the same fraction as in the universe 
with a bias factor b. We substitute this assumption in 
Eq. (66) to achieve 



fgas 



(68) 



Using the distribution of gas and matter in cluster, Sasaki 
(1996) showed that the fraction of gas depends on angu- 
lar distance with f gas oc D^ 2 [20]. On the other hand 
the fraction of gas obtained from the observation depends 
on the model we are assuming for the dynamics of the 
universe. It is assumed that f gas should in reality be in- 
dependent of the rcdshift. To determine the constraints 
on the proposed modified gravity action, we fit the f gas 
data with a model that accounts for the expected ap- 
parent variation in f gas {z) as the underlying cosmology 
is varied. We choose both SCDM (a flat universe with 
TO = 1 and h = 0.5) and ACDM as reference cosmology 
models. The ratio of gas fraction for a given model to the 



reference model is: fgas^/fg 
On the other hand using Eq. (68) for a given model, the 
gas fraction for a reference model is obtained by: 



/{mod) 
gas 



D 



(ref) 



ID 



(mod) 
A 



3/2 



f(ref) _ 
J Q 



^ 



I gas 



3/2 



(69) 



where superscrips (ref) correspond once to SCDM and 
then to the ACDM model [25]. We use x 2 to compare 
gas fractions of observational and theoretical models as 
follows: 
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FIG. 11. Marginalized likelihood functions of f(R) modi- 
fied gravity model free parameter, n. The solid line corre- 
sponds to the joint analysis of SNIa data (new gold sam- 
ple) and CMB, the dashed line shows the joint analysis 
of SNIa+CMB+SDSS data,the dash-dot line corresponds 
to SNIa+CMB+SDSS+LSS and dash-dot-dot line indicates 
SNIa+CMB+SDSS+LSS+LCDM. The intersections of the 
curves with the horizontal solid and dashed lines give the 
bounds with la and 2a level of confidence respectively. The 
results for two former analyses are very similar. 



from the WMAP, recently observed baryonic peak from 
the SDSS and 2dF and the gas mass fraction in cluster 
of galaxies to constrain the parameter of the modified 
gravity model by minimizing the combined \ 2 = XsNia + 



Xcmb + Xbao 



Y 2 

Agas 



The best values of the parameters of the model from 
the fitting with data, including SNIa new sample are 

n = 0.9118:8?, X = 3.67ig:« and ^« = 0.29±g;SB 
and data with including SNIa union sample results in 
n = 0.9818:81, X = 4.39l°;l and a m = 0.25l°;J? - 
Here we marginalized overall Hubble parameter in likeli- 
hood analysis. Figures (11) and (12) show the likelihood 
function as function of exponent n. Also Figures (13) 
and (14) represent the likelihood function of X in two 
different supernova data sets. 



VI. CONSTRAINTS BY LARGE-SCALE 
STRUCTURES:DYNAMICAL PARAMETER 



Xgas 



( fobs fthe\2 

\J qas J qas ) 



gas 



This is the fourth geometrical constraint. 



(70) 



E. Combined analysis: 
SNIa+CMB+BAO+GAS-FRACTION 

In this section we combine SNIa data (from SNIa new 
Gold sample and mixed SNLS), CMB shift parameter 



So far we have only considered the observational re- 
sults related to the background evolution. In this section, 
using the linear approximation of structure formation, 
we obtain the growth index of structures and compare 
it with the result of observations by the 2-degree Field 
Galaxy Rcdshift Survey (2dFGRS). As we mentioned be- 
fore, the evolution of the structures depends on both the 
dynamics of the background and the gravity law that 
governs the dynamics of particles inside the structure. 

Here the evolution of structures in the modified grav- 
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FIG. 12. Marginalized likelihood functions of f(R) 
modified gravity model free parameter, n. The solid 
line corresponds to the joint analysis of SNIa data 
and CMB, the dashed line shows the joint analysis of 
SNIa+CMB+SDSS data, dash-dot line corresponds to 
SNIa+CMB+SDSS+LSS and the dash-dot-dot line indicates 
SNIa+CMB+SDSS+LSS+LCDM. The intersections of the 
curves with the horizontal solid and dashed lines give the 
bounds with la and 2a level of confidence, respectively. The 
results for two former analysis are very similar. 
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FIG. 13. Marginalized likelihood functions of f(R) mod- 
ified gravity model as a function of Xq. The solid line 
corresponds to the joint analysis of SNIa data (new Gold 
sample) and CMB, the dashed line shows the joint analy- 
sis of SNIa+CMB+SDSS data, dash-dotted line corresponds 
to SNIa+CMB+SDSS+LSS and dash-dot-dot line indicates 
SNIa+CMB+SDSS+LSS+LCDM. The intersections of the 
curves with the horizontal solid and dashed lines give the 
bounds with la and 2a level of confidence respectively. The 
results for two former analysis are very similar. 



ity will be studied through the spherical collapse model. 
Recently a procedure has been put forward by Lue, Scoc- 
cimarro, and Starkman (2004) which relies on the as- 
sumption that Birkhoffs theorem, holds in a more gen- 
eral setting of modified gravity theories. This procedure 
also is applied in the Palatini formalism of f(R) gravity 
[26]. According to this procedure, it is assumed that the 
growth of large-scale structure can be modeled in terms of 
a uniform sphere of dust of constant mass. This structure 
evolves as a FRW universe. Using Birkhoffs theorem, the 
space-time at the empty exterior of this structure is then 
taken to be a Schwarzschild-like metric. The components 
of the exterior metric are then uniquely determined by 
smoothly matching the interior and exterior regions. In 
the Palatini formalism the metric outside the spherical 
distribution of matter depends on the density of matter 
which may modify the Newtonian limit of these theories, 
however here we assume Schwarzschild-like Newtonian 
limit. 

The continuity and Poisson equations for the density 
contrast 8 = 8pj p in the cosmic fluid provides the evolu- 
tion of density contrast in the linear approximation (i.e. 
8 < 1) as 



8 + 2-8- AirGp 8 
a 



0, 



(71) 



where the dot denotes the time derivative and we assume 
the size of structures to be larger than the Jeans length. 
The effect of a modified gravity results from the modi- 
fication to the background dynamics and we adopt the 



same Poisson equation for the weak field regime. In or- 
der to use the constraint from large-scale structure, we 
rewrite the above equation in terms of X. So we have 



d 2 S 
da^ 



dS_ 

da 



n'{x) dx 

H{X) da 



2H 2 (X)w 



■8 = (72) 



In the standard linear perturbation theory, the peculiar 
velocity field v is determined by the density contrast as 



v(x) = H Q 



J_ 

4tt 



rd 3 y, 



where the growth index / is defined by 



dhxS 
d In a 



(73) 



(74) 



and it is proportional to the ratio of the second term of 
Eq. (71) (friction) to the third term (Poisson). 

We use the evolution of the density contrast 5 to com- 
pute the growth index of structure /, which is an im- 
portant quantity for the interpretation of peculiar veloc- 
ities of galaxies. Replacing the density contrast with the 
growth index in Eq.(73) results in the evolution of growth 
index as 



3a 



df _ 
din a 2aH 2 (X) 



-f-f 



aH'(X) dX 



H(X) da 



(75) 



To put constraint on the model using large structure 
data, we rely on the observation of 220000 galaxies with 
the 2dFGRS experiment, which provides a numerical 
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FIG. 14. Marginalized likelihood functions of f(R) modi- 
fied gravity model as a function of Xo. The solid line corre- 
sponds to the analysis of SNIa data and CMB, the dashed line 
shows the joint analysis of SNIa+CMB+SDSS data, dash-dot- 
ted line corresponds to SNIa+CMB+SDSS+LSS and dash- 
dot-dot line indicates SNIa+CMB+SDSS+LSS+LCDM. The 
intersections of the curves with the horizontal solid and 
dashed lines give the bounds with la and 2a level of con- 
fidence respectively. The results for two former analysis are 
very similar. 

value for the growth index. By measurements of the two- 
point correlation function, the 2dFGRS team reported 
the redshift distortion parameter of (3 = f jb = 0.49±0.09 
at z = 0.15, where b is the bias parameter describing the 
difference in the distribution of galaxies and their masses. 
Verde et al. (2003) used the bispectrum of 2dFGRS 
galaxies [27,28] and obtained b ver d e = 1.04 ± 0.11 which 
gave / = 0.51 ± 0.10. Now we fit the growth index at 
z = 0.15 derived from the Eq.(75) with the observed 
value. 

2 [fobs(z = 0.15) -fth(z = 0.15; X )} 2 
Xlss = -2 ■ ( 76 ) 

a fobs 

Finally we do likelihood analysis with considering all 
the observations and obtain the 2D distribution of a like- 
lihood function in terms of n and Xq in Fig. (17). 



A. Perturbation Theory 

In the previous section we have seen the effect of mod- 
ified gravity on the structure formation in the weak filed 
regime through the background effect. In this method 
changing the dynamics of universe (i.e. scale factor) al- 
ters the formation of the large-scale structures. 

In this section we study the relativist structure forma- 
tion theory through the perturbation in the homogenous 
background metric and energy-momentum tensor. This 
approach may uncover whether modified gravity theo- 
ries driving late-time acceleration predict any testable 



features on CMB or large-scale structures in linear or 
nonlinear regimes. 

Let us consider a flat universe dominated by pressure- 
less cold dark matter. We identify perturbation in con- 
formally flat FRW space-time by ten elements as follows: 



a 2 (^){-(! + 2a)dr/ 2 - 2(/3 l + b^drjdx 1 



+[$ + %gV + %ij + c m ) + %]d*W} (77) 

where i] is the conformal time. cv„3,7,and <f> are the scalar 
perturbations to the metric, bi and Ci are divergenceless 
vectors where each one with 2 degrees of freedom and Ziy- 
is a traceless-divergenceless symmetric 3x3 matrix with 
2 degree of freedom. 

On the other hand perturbation of conservation of 
energy-momentum tensor results in the continuity and 
Euler equations as follows: 



5 = —kv + an — 3Ha, 
v = —Hv + fca, 



(78) 
(79) 



where k — 3a 1 (Ha 



l Ox and x = a (P + 7). 



Using conformal gauge, the field equation for the den- 
sity contrast in the modified gravity framework obtain as 
[29]: 



(80) 



where ' is derivative with respect to the conformal time 
and H — — = a. f and C arc defined as 

a ^ ^ 



f = 1 

C = i 



2FF H~2F H- 2FF H 
FH 2 (2FH + F') 

(i-O 



H" - 2H'H 
FH 



3(2FH + F')(H" -2H'H) 



(81) 



(82) 



where k is the wave number of structures in the universe. 
In the case of the Einstein-Hilbert action, F' = F" = 
which results in £ = 1 and subsequently £ = 1. The dif- 
ferential equation governing the evolution of the density 
contrast in this case reduces to: 



5" + HS' - (— - 2H)S = 0. 
H 



(83) 



For comparison of Eqs.(80)and (83), we obtain the differ- 
ence in the density contrast between the ACDM model 
and the modified gravity as indicated in Fig. (15) for a 
structure with the size of k = O.OlMpc -1 . For the larger 
scales, the third term in the Eq.(82) tends to zero and 
we get smaller difference between the density contrast 
in these two solutions. In order to compare these re- 
sults with data, more detailed simulation in the nonlinear 
regime of the structure formation is essential. 




FIG. 15. Difference between the density contrast in ACDM 
model and Modified gravity model. For larger redshifts the 
difference between these two solutions is negligible. 



VII. AGE OF UNIVERSE 



from the combined analysis of SNIa+CMB+SDSS+LSS 
is 14.69±g || Gyr and 13-45±g;|g Gyr for new gold sample 
and union data sample, respectively. These values are in 
agreement with the age of old stars [30]. Here we take 
three OHRG for comparison with the modified gravity 
model considering the best fit parameters, namely the 
LBDS 53W091, a 3.5-Gyr old radio galaxy at z = 1.55 
[33], the LBDS 53W069 a 4.0-Gyr old radio galaxy at 
z = 1.43 [34] and a quasar, APM 08279+5255 at z = 3.91 
with an age of t = 2.1^° iGyr [35]. To quantify the age- 
consistency test we introduce the expression r as: 



t(z;X ) = t(z;X )H Q 

tobs t b s Ho 



(85) 



where t(z) is the age of universe, obtained from the 
Eq.(84) and t b s is an estimation for the age of an old 
cosmological object. In order to have a compatible age 
for the universe we should have r > 1. Table III reports 
the value of r for three mentioned OHRGs with various 
observations. We see that f(R) modified gravity with the 
parameters from the combined observations, provides a 
compatible age for the universe, compared to the age 
of old objects, while the SNLS data result in a shorter 
age for the universe. Once again, APM 08279 + 5255 at 
z = 3.91 has a longer age than the universe but gives 
better results than some of modified gravity models [36] . 



The age of the universe integrated from the big bang 
up to now for a flat universe in terms of free parameters 
of model n and Xq is given by: 



to{X p ) 



dt = 



dz 



3Ho 



o 
F' 



(l + z)H(z) 
- XF" dX 



X,; 



2F - XF' H(X) 



(84) 



Figure (16) shows the dependence of H t (Hubble pa- 
rameter times the age of universe) on Xq for a flat uni- 
verse. In the lower panel we show the same function for 
ACDM universe in terms of f^A for comparison. As we 
expected, Xq in modified gravity behaves as a dark en- 
ergy and increasing it makes a longer age for the universe, 
in the same direction as increasing the cosmological con- 
stant. 

The "age crisis" is one the main reasons for the accel- 
eration phase of the universe. The problem is that the 
universe's age in the CDM universe is less than the age 
of old stars in it. Studies on the old stars [30] suggest an 
age of 13^2 Gyr for the universe. Richer et. al. [31] and 
Hasen et. al. [32] also proposed an age of 12.7 ±0.7 Gyr, 
using the white dwarf cooling sequence method. 

We use the age of universe in this model for the con- 
sistency test and compare the age of universe with the 
age of old stars and old high redshift galaxies (OHRG) in 
various redshifts. Table II shows that the age of universe 



VIII. CONCLUSION 

In this work we proposed the action of f(R) = (R n — 
Rq ) 1 /" to obtain the dynamics of the universe. We used 
the Palatini formalism to extract the field equation. The 
advantage of this formalism is that the field equation is 
a second-order differential equation and in the solar sys- 
tem scales we can recover a Schwarzschild-de-Sitter space 
with an effective cosmological constant compatible with 
the observations. The other advantage of the Palatini 
formalism is that it does not suffer from the curvature 
instability as pointed out in [37]. 

We used cosmological tests based on background dy- 
namics such as Supernova Type la, CMB-shift parame- 
ter, baryonic acoustic oscillation and gas mass fraction 
of the cluster of galaxies. We also used data from the 
structure formation to put constrains on the parameters 
of the model. Table II represents constrains on the pa- 
rameters of model considering the observational data and 
their combination. We also showed that this model pro- 
vides an age for the universe sufficiently longer than the 
age of old astrophysical objects. 

Comparing this model with the observations we put 
the constrain of n = 0.98lo os f° r ^ ne exponent of action 
and X = 4.39±£ij| or equivalently Cl m = 0.25lg;^. The 
best value for this model shows that a standard ACDM 
model also reside in this range of solution. Our result 
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FIG. 16. Hoto (age of universe times the Hubble constant 
at the present time) as a function of Xq (upper panel). H t 
for ACDM versus Q\ (lower panel). Increasing Xo gives a 
longer age for the universe. This behavior is the same as 
direction as in ACDM universe. 
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FIG. 17. Joint likelihood function in terms of n and Xo 
considering all the observable data. In the upper panel the 
SNIa data is taken from gold sample and in the lower panel 
the data is Union sample. 



is in agreement with the recent work by Kowalski et al. 
(2008) where they also obtained almost a ACDM universe 
with a nearly constant equation of state for a dark energy 
model [38]. 
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TABLE II. The best values for the parameters of modified 
gravity with the corresponding age for the universe from fit- 
ting with SNIa from new Gold sample and Union data sample, 
SNIa+CMB, SNIa+CMB+BAO , SNIa+CMB+BAO+LSS 
and SNIa+CMB+BAO+LSS+GAS(ACDM) experiments at 



one and two a confidence level. The value of Q m is deter- 
mined according to equation (25) 
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TABLE III. The value of r for three high redshift objects, 
using the parameters of the model derived from fitting with 
the observations at one and two a level of confidences. 
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